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Abstract. We present our view of the main physical ingredients determining the evolution of 
neutron star magnetic fields. This includes the basic properties of neutron star matter, possible 
scenarios for the origin of the magnetic field, constraints and mechanisms for its evolution, and 
a discussion of our recent work on the Hall drift. 


INTRODUCTION 


In the present eonferenee, Mariano Mendez reviewed our observational knowledge 
about neutron star magnetie fields, while Don Melrose diseussed eurrent models of the 
pulsar magnetosphere. Therefore, our review foeuses on theoretieal ideas regarding the 
magnetie field inside neutron stars, paying partieular attention to proeesses by whieh 


this field might evolve in time. We refer to our previous reviews II28L13011 for more basie 
and general diseussions of our knowledge of neutron star magnetie fields. 

Our diseussion starts with an overview of the properties of matter in neutron star, 
whieh are obviously an essential ingredient in any model for the evolution of the mag¬ 
netie field. We eontinue with a brief diseussion of the state of knowledge regarding the 
origin of the magnetie field from the main-sequenee stage through the hot, protoneu¬ 
tron star phase, whieh sets the initial eonditions for the subsequent evolution. Then, 
we quiekly present the observational evidenee for the evolution of the magnetie field, 
followed by the basie equations that are likely to deseribe this evolution and the main 
proeesses involved in it. One of these proeesses, the Hall effeet, is singled out for more 
detailed diseussion, where we present reeent, yet unpublished analytie solutions for sim¬ 
ple geometries, and diseuss their implieations. 



MATTER IN NEUTRON STARS 


Protoneutron stars are bom through gravitational eore eollapse in massive stars, whieh 
leaves them at a high temperature, T ~ lO^^K lOMeV, eomparable to the Fermi 
energies of their main eonstituent partieles. However, neutrino emission ensures rapid 
eooling, so all eurrently observed neutron stars (with ages lO^yr or higher) are 
expeeted to have internal temperatures at least two orders of magnitude lower, at whieh 
the matter is highly degenerate and likely elose to its quantum ground state, although 
small perturbations from the latter may eause interesting effeets ll^ [pl [T11l . some of 
whieh are important for the evolution of the magnetie field Eilllllli and will be 
diseussed in what follows. 

In this highly degenerate, near-equilibrium state, the properties of matter are deter¬ 
mined by a single parameter, for example density or pressure. The density inside neutron 
stars eovers many orders of magnitude, allowing for different regimes in the properties 
of matter, as diseussed below. For a mueh more detailed diseussion, see Ref. 11^ . 


The crust 


The neutron star erust is, by definition, the outer part of the star, where atomie nuelei 
are present. As in ordinary matter, they eoexist with eleetrons, whieh, due to the high 
density, are too energetie to be bound to individual nuelei. At T ~ lO^^K, the nuelei are 
expeeted to freeze into a solid, leaving the eleetrons as the only moving eharge earners. 

Traditionally, it has been assumed that the erust is in its absolute ground state, in 
whieh one partieular kind of nueleus, whieh minimizes the enthalpy at the loeal pressure, 
is present at any given depth in the star, and these nuelei form a near-perfeet erystal 
lattiee. This pieture has been ehallenged by Jones lllSh . who argued that thermodynamie 
fluetuations at the time of freezing would ensure the presenee of several kinds of nuelei 
at any given pressure, yielding a very impure solid. This result may be very important 
for the evolution of the magnetie field, as it would substantially inerease the resistivity 
of the neutron star erust, partieularly at low temperatures M. 

^gem 

the erust. At a temperature similar to that for freezing the lattiee, these neutrons are 
expeeted to form Cooper pairs and turn superfluid, largely deeoupling dynamieally from 
the rest of the star. This has been invoked to explain pulsar “glitehes”, quiek inereases of 
the rotation rate of the externally visible parts of a neutron star, whieh might be eaused by 
a eatastrophie transfer of angular momentum from the more rapidly rotating superfluid 
neutrons. Given that these neutrons are uneharged and internet only weakly with the 
eleetrons, their presenee is not likely to be relevant for the magnetie field evolution. 


At densities above “neutron drip”, ~ 4 x lO^^gem free neutrons also appear within 


The core 

At ~ 2 X 10^‘^gem slightly below standard nuelear density, the nuelei are expeeted 
to loose their individual identity, so matter at higher densities eonsists of a liquid-like 




mix of neutrons (n), protons (p), and electrons (e^), which at progressively higher 
densities are joined by additional particles, such as muons, hyperons, and mesons. 
Chemical equilibrium among all these particles is established by weak interactions such 
as the neutron beta decay (n ^ p + e^ + v) and inverse beta decay (p + e^ ^ n + v), 
in most cases involving the emission of neutrinos (v) or antineutrinos (v), which easily 
leave the star. In this phase, the strong Pauli blocking of final states reduces the cross 
sections for, e.g., electron-proton collisions, yielding a very high electrical conductivity 
iQ]. Here, all particles are movable (within constraints to be discussed below), making 
the evolution of the magnetic field much more complicated. 

Models predict that the strongly interacting particles in this phase (n, p, possibly hy¬ 
perons) again form Cooper pairs and enter a superfluid state, in which the neutron vor- 
ticity is concentrated in quantized vortex lines of microscopic thickness, much smaller 
than their average spacing, and the magnetic flux may be similarly concentrated in pro¬ 
ton vortices. The transition temperatures for these superfluid states are highly uncertain. 
Their effect on the thermal evolution of the neutron star, although potentially quite im¬ 
portant, is not obviously shown by the observations of cooling, young neutron stars (e.g.. 
Fig. 1 of Ref. ll^ l. although better fits to the (quite uncertain) data have been obtained 
by allowing for neutron stars with a range of superfluid properties (with mass as the 
controlling, free parameter) ll^ 1^ . On the other hand, the core superfluids are not 
expected to be as easily decoupled as the crustal neutron superfluid, and therefore play 
no role in glitches. Therefore, we consider it safe to say that there is no strong evidence 
for the existence of superfluids in the neutron star core. This conclusion, in addition to 
the fact that superfluids introduce a new element of uncertainty in the already highly 
complicated models of magnetic field evolution, motivate us to ignore it in most of what 
follows, although it has been discussed in detail by other authors 


ORIGIN 

The range of observed surface magnetic field strengths in degenerate stars is quite large, 
from below lO'^G to fimax ~ lO^G on white dwarfs, and from 10^ ^G up to fimax ~ lO^^G 
for the dipole fields on young neutron stars, including classical radio pulsars, anomalous 
X-ray pulsars, and soft gamma-ray repeaters (the field on millisecond pulsars may have 
been decreased by processes related to accretion, and not be related to the original field). 
The measurable dynamic range on massive main sequence stars is much smaller, but it 
is suggestive that the largest inferred magnetic fluxes in all three kinds of objects are 
similar, $ = nR^Bmax ~ 3 x lO^^Gcm^, where R is the stellar radius . This implies 
that, in principle, flux freezing from the main sequence is enough, and no additional 
magnetic field generation mechanisms need to be invoked to explain even the strongest 
fields observed on the compact stars. 

On the other hand, stellar evolution, specially in the massive stars that eventually 
become neutron stars, is quite eventful, so the field geometry (and perhaps its strength) 
is likely to become modified along this evolution, possibly coupled to the evolution of 
the angular momentum. In particular, the vigorous convection occurring in a newborn, 
perhaps rapidly rotating, protoneutron star is an ideal setting for a dynamo to operate. 


which might increase even a fairly weak initial field up to IQi^GlUl. 

On the other handj thermomagnetie instabilities driven by the heat flow through the 
neutron star erust 11391 14 l^h ean at most account for pulsar-like field strengths IQi^G, 
making them mueh less likely to be relevant. 

Thus, although no detailed seenario has been eonvineingly worked out, it appears 
likely that a strong magnetie field, permeating mueh of the core of the star, is present at 
least from the hot, protoneutron star phase onwards. 

Regarding its geometrieal eonfiguration, reeent magnetohydrodynamie (MHD) simu¬ 
lations by Braithwaite & Spruit [j] show that, in a stably stratified star, a eomplieated, 
random, initial field generally evolves on an Alfven-like time scale to a relatively sim¬ 
ple, roughly axisymmetric, large-seale configuration containing a toroidal and a poloidal 
component of comparable strength, both of whieh are required in order to stabilize eaeh 
other. Only the poloidal eomponent erosses the surfaee of the star, yielding an exter¬ 
nal field that is roughly dipolar, though perhaps somewhat offset from the eenter, as is 
observed in magnetic A stars and white dwarfs. 

We should note, however, that both A stars and white dwarfs are of essentially uni¬ 
form eomposition (exeept for the diseontinuity between the hydrogen-burning, eonvee- 
tive eore and the ehemieally pristine radiative envelope in A stars), so they are stably 
stratified only through the entropy gradient in their radiative regions. This means that, 
if stable stratifieation is an essential ingredient in stabilizing the configurations found 
in Ref. |0], their lifetime would be limited by the exehange of entropy between regions 
of different field strengths O. If, as seems to be the ease, this time is shorter than 
the lifetime of the stars, the field eould be eonfined only by eompositional gradients 
(discussed below for the case of neutron stars), as in the core-crust transition of upper 
main-sequence stars. 

In any ease, it appears plausible that fields of eonfigurations similar to those of Ref. 
10], with poloidal field strengths ~ 10^^~^^G, and perhaps with their toroidal eomponent 
enhanced by the initial differential rotation 113811 . could be present in early neutron stars. 


EVIDENCE FOR MAGNETIC FIELD EVOLUTION 

Sinee the earliest days of the study of pulsars, it has been claimed JUl lISl that the 
distribution of these objeets on the period (P) - period derivative (P) plane (sometimes 
eombined with other variables, sueh as their position with respeet to the Galaetie plane, 
their proper motion, and their luminosity) would indieate a magnetie field deeay on a 
time scale eomparable to a pulsar lifetime. A safe proof or disproof of such statements is 
very diffieult to make, sinee various seleetion effeets are important, but recent work has 
generally not eonfirmed them |0,l3l ■ Nevertheless, in the past it provided the motivation 
for a fair amount of theoretieal work trying to explain the elaimed deeay. 

Another weak argument for magnetic field evolution is the fact that the braking indiees 
n = measured in very young pulsars have always turned out to be smaller than 

the value (n = 3) expected from a magnetie dipole torque ll^ with a eonstant dipole 
moment. This may indieate that the magnetie field inereases in very young pulsars, or 
else that the dipole formula is not eompletely adequate (and therefore does not give a 



very reliable estimate of the dipole field strength lEl). 

A stronger, “elassieal” argument for field deeay is the observation that old neutron 
stars, such as found in low-mass X-ray binaries and millisecond pulsars, tend to have 
much weaker fields (~ 10^“^G) than the younger classical pulsars and high-mass X- 
ray binaries (~ lO^^G). It is still not clear if this is spontaneous decay happening 
over the long lifetime of these objects^or whether it is due to the accretion process 
diamagnetically screening the field 11331 19|] or increasing the resistivity to cause a faster 
decay. In all these cases, it seems surprising that the field in all millisecond pulsars 
appears to reach a bottom value ~ lO^G and does not decay beyond that. A possible 
explanation might be that the original field of the star is completely screened or lost 
from the star, whereas new magnetic flux is carried onto the neutron star by the accretion 
flow, limited by the condition that the fluid stresses be strong enough to compress the 
field onto the star, which yields a maximum field strength of about the right order of 
magnitude. 

Finally, but most interestingly, soft gamma-ray repeaters and anomalous X-ray pulsars 
appear to radiate substantially more power than available from their rotational energy 
loss. At the same time, their dipole field inferred from the observed torque is stronger 
than in all other neutron stars, ~ lO^^^^^G (although a couple of pulsars with B ~ lO^^G 
are also known illzillii). This makes it natural to accept the argument of Thompson 
& Duncan 10 that these sources are in fact “magnetars”, powered by the dissipation 
of their magnetic energy (which requires an rms internal field even stronger, though not 
by very much, than the inferred dipole field). In addition, torque changes of both signs, 
associated with outbursts from these sources, also argue for a changing magnetic field 
structure. These are, at the moment, the only objects in which there is strong evidence 
for spontaneous evolution of the magnetic field. 


A PHYSICAL MODEL 

As explained above, a neutron star core contains a number of mobile particle species. 
After an initial transient, following the formation of the neutron star, on which all sound 
and Alfven waves are damped, the evolution of the magnetic field should be slow enough 
to make the inertia of the particles negligible, leading to the equation of diffusive motion 
for the particles of each species i, 

0 =x.B^ -Y,7iinj{vi-Vj), (1) 

where jii, m*, qu Vi are their chemical potential (Fermi energy), effective mass (including 
relativistic corrections due to random motions and interactions), electric charge, and 
mean velocity, i// is the gravitational potential, E and B are the electric and magnetic 
fields, and the last term represents the momentum transfer due to collisions with all 
other particle species j, each with number density nj. The collisional coupling strengths 
are parameterized by the symmetric matrix jij, whose coefficients generally depend on 
position. 


The collision terms have two effects in the context of the evolution of the magnetic 
field: 

(a) They damp the relative motion of positive and negative charge carriers, leading to a 
resistive diffusion of the magnetic field, opposed by the induced electric field. For a 
magnetic field of spatial scales comparable to the neutron star radius, the induction 
is large and the resistance is small, as in most astrophysical settings, so this process 
is slow and unlikely to have any effect over the lifetime of any neutron star It 
may, however, be important if the magnetic field is created within a thin surface 
layer (which does not fit with the likely formation scenarios we described) or if 
small-scale structure is created by other processes, as discussed below. 

(b) Collisions tend to keep the different particle species moving with similar velocities, 
as in the standard MHD approximation. We now proceed to analyze this type of 
motion. 

Adding the forces given by eq. dH) over all particles (of all species) within a volume 
containing a unit total mass and no net charge, one obtains the equation of MHD 
equilibrium. 


^ VF ^ jxB 

0 =-, 

p pc 


( 2 ) 


where the mass density p = the current density j = = (c/4;r)V x B, 

and the pressure gradient term was obtained from the zero-temperature Gibbs-Duhem 
relation, dP = Taking the curl of eq. Q, we obtain 


vfxvp^v 

P \ / 

The right-hand side is generally non-zero, therefore equilibrium is only possible if 
the pressure and density gradients appearing on the left-hand side are not parallel. 
This is not possible in cold matter in chemical equilibrium, since then the density is 
a unique function of pressure (because all the chemical abundances are also determined 
by the latter). Therefore, the magnetic field generally perturbs the chemical equilibrium, 
causing a tiny misalignment between the pressure and density gradients. 

A magnetic force distribution that has a horizontal curl component induces compen¬ 
sated up and downward motions of the fluid in different regions, producing the men¬ 
tioned misalignment of the density and pressure gradients, which chokes the motion. 
Only the horizontal fluid motions produced by a magnetic force density with a purely 
vertical curl can proceed essentially unimpeded. This is a manifestation of the stable 
stratification of the neutron star matter 11^ 1^ 12911 . which prevents even a magnetar- 
strength field from being transported vertically through the star. The only ways by which 
this constraint can be circumvented are 11241 iT^ : 

(a) elimination of the induced chemical imbalance by weak interaction processes, 
which is most effective at high temperatures, where reactions are fast, and 

(b) relative motion of different particle species. 





We now proceed to discuss each of these processes. 


Bulk motion facilitated by weak interactions 

As described above, the Lorentz force produces non-barotropic perturbations to the 
pressure, of magnitude 8P Y, Hid Hi, where the relation to the chemical 

potential perturbations again comes from the Gibbs-Duhem relation. If, for definiteness, 
we consider the simplest possible neutron star matter, composed of neutrons, protons, 
and electrons (the latter two related by the condition of charge neutrality, Up = = ric), 

the total chemical imbalance is 


AjU = \iip + He-Bn\ r^B^/{%Ttnc) ~ Sfifj keV. (4) 

where 5i5 = 5/(10^^G). An asymmetry in the weak interaction rates tends to reduce 
this imbalance. In what follows, we assume that the dominant process are “modified 
Urea reactions” without Cooper pairing, which give a reasonable fit to the early cool¬ 
ing of neutron stars (see Fig. 1 in Ref. lH^]) and possibly to the late, “rotochemical” 
reheating of millisecond pulsars im. As long as r 3> A/i ll^ . the imbalance decays 
exponentially, with time constant tmU ~ 0.5/T^ yr, essentially the cooling time of the 
star, which is highly sensitive to temperature {T = Tg x lO'^K^TgxSe keV). In this 
time, the Lorentz force makes the fluid move a small fraction, ~ A/r //tg ~ 2 x 10 - 5523 , 
of the stellar radius, so the total decay time of the field is 

^l, 3 X 10* 

yr, (5) 

in principle quite short for strong fields and high temperatures. However, it is important 
to note that it is longer than the cooling time by the generally quite considerable ratio 
Thus, unless the magnetic field is in the vicinity of ~ 10 ^^ q (making A/i ~ He)-, 
a heat source is required to prevent passive cooling of the neutron star from “freezing” 
the magnetic field. The most obvious energy source is the magnetic field itself yj], 
whose energy in this scenario is released by the weak interactions. We note, however, 
that, in the absence of a chemical imbalance, the neutrino emission produced by these 
weak interaction leads to the cooling of the star. It is balanced by the heat release only if 

the chemical imbalance is fairly large, A/i k, 5.5T lEl] or, equivalently, 5i5 135^'^^. 

A highly magnetized neutron star would be born at a high temperature and quickly cool 
down until this condition is satisfied, after which the temperature would be stabilized 
(at Tg ~ 0.2[10^ yr/t]^/2) by the energy injected through the decaying magnetic field,' 
5 i 5 ~ 6.5(10^ yr/t) At this fairly high temperature, the collision rates are high, thus 
the single-fluid MHD approximation is probably adequate. 


' The coefficient is slightly higher than found by Thompson & Duncan It^ . who first obtained essentially 
the same result. 




Relative motion of charged and neutral particles 


The process just described is not relevant for stars with internal fields substantially 
lower than 5 x 10^^ G, which would enter the photon-cooling epoch before the 
magnetic field has substantially decayed or contributed to any reheating. For these, 
the only possibility to have magnetic field decay relies on relative motion of charged 
particles. We consider a slight extension of the model of npe matter proposed by 
Goldreich & Reisenegger llT^ . in which different linear combinations of the eqs. o 
for the three particle species, together with the induction equation dB/dt = —cV x E, 
yield 



xV(^p + jUg). (6) 


In the absence of the last two (resistive and battery) terms, which in realistic conditions 
are generally small, the magnetic field is advected by a combination of three velocities: 

(a) the neutron velocity v„, corresponding to the bulk motions considered in the previ¬ 
ous subsection; 

(b) the ambipolar diffusion velocity. 


YpniVp Vfi)Yen{Ve Vn) _ j^B/(^n^c') V(A^) 

Ypn~Y Yen T U^;) -(- Yen') 


(7) 


va = 


representing a relative motion of the charged particles with respect to the neutrons, 
driven by the Lorentz force, and potentially choked by the gradient it induces in the 
chemical potential imbalance; and 

(c) the Hall drift velocity. 



( 8 ) 


due to the relative motion of protons and electrons and proportional to the current 
density. 

The chemical imbalance induced by ambipolar diffusion can in principle be elimi¬ 
nated by weak interactions, on time scales similar to those of the bulk motions, which 
are short only at high temperatures, at which ambipolar diffusion is strongly suppressed 
by collisions. The regime in which it may be important is at lower temperatures, where 
collisions are less constraining, but weak interactions are strongly suppressed. In this 
case, only the solenoidal mode of ambipolar diffusion (driven by the finite-curl, zero- 
divergence part of the Lorentz force 111 211 1 can proceed (and even this only in the sim¬ 
plified case of a uniform charged fluid with Up = Ug, without the stabilizing effect of 
nonuniform relative particle abundances if more charge carriers are present). 

The Hall drift is of a very different character, since the drift velocity is directly related 
to the magnetic field, and independent of any arising fluid forces. In a solid medium 
such as the neutron star crust, where the electrons are the only free charge carrier, it 









is the only active process besides resistive diffusion. This regime has been studied by 
several authors and is discussed in the following section. The more complicated regime 
of the fluid core, in which the Hall effect acts in conjunction with (or in opposition to) a 
highly constrained ambipolar diffusion, is largely an open question [0]. 


HALL DRIFT 


In the neutron star crust or another solid, conducting medium, the only moving charges 
are the electrons, and the evolution of the magnetic field is described by the “Hall 
equation”. 


dt 


Vx 


( - - —V X 5^ X 5+ -^V X B 

\ Annee ) Ako 


(9) 


where the two terms on the right-hand side correspond to the Hall drift and resistive 
diffusion, respectively, whose relative importance in the neutron star crust is still a matter 
of controversy llKlIl^ . 

Here, we discuss some physical issues that are likely to determine the evolution of 
magnetic field under the Hall effect, without attempting to cover the many simulations 
recently performed to address how specific magnetic field configurations might decay in 
real neutron stars SEi. _ 

Goldreich & Reisenegger focused on the more interesting case in which, on large 
scales (comparable to the crust thickness) the Hall effect is dominant, and argued, by 
analogy with the Euler equation of fluid dynamics, that the nonlinear Hall term may 
give rise to a turbulent cascade to small scales. Due to the presence of stable, linear 
modes, this turbulence would be “weak”, with an energy transfer time generally longer 
than the typical oscillation period on a given scale, resulting in a power spectrum oc 
slightly different from the Kolmogorov spectrum oc of fluid turbulence. At small 
scales, the magnetic energy would finally be dissipated by resistivity. This nonlinear 
evolution was simulated by Biskamp et al. |S|, who found an energy cascade to small 
spatial scales, but with an even steeper spectrum than predicted, 

A complementary approach is to study analytic solutions of the Hall equation (with or 
without the resistive term), which has been done by Vainshtein et al. Il40ll . Gumming et 
al. and by our group HQ . In the latter, which extends and generalizes the work 
of Ref. 1^ . we first considered a purely toroidal field, B = where R,^,z 

are the standard cylindrical coordinates. When evolved by the Hall equation, the field 
remains toroidal. In order to describe its evolution, it becomes convenient to introduce 
a new coordinate x = c/[Anene{R,z)R^]. Surfaces of constant x (hereafter, ;(^-surfaces) 
are toroids contained inside the star. A complementary coordinate s can be defined on 
each ;f;-surface by d/ds = —R^V^ x Vx ' terms of which the Hall equation (now 
neglecting the resistive term) reduces to the Burgers equation. 






ds 


0 , 


( 10 ) 


with the well-known, implicit solution = f{s — ^t). Each value of ^ is carried 
around the corresponding ;t^-surface with a velocity proportional to developing 








discontinuities at the (eomoving) points where jds is large. There, the resistive term 
inereases and tends to smooth the diseontinuity. Therefore, magnetie energy is dissipated 
at the (relatively fast) rate as it is fed into the discontinuity by the Hall drift. Eventually, 
on the eharaeteristie time seale of the Hall drift, beeomes uniform on eaeh ;K-surfaoe, 
so the resulting magnetic field B = evolves on the (assumed) much longer 

resistive time seale. 

However, further analysis ll^ shows that this field is unstable to small, poloidal per¬ 
turbations, which grow when different segments of the poloidal field lines, crossing dif¬ 
ferent ;t^-surfaces, are earried along at different speeds by the Hall drift veloeity associ¬ 
ated to the toroidal field. So, the question arises: Are there any non-trivial eonfigurations 
whieh are not only statie, but also stable, under the Hall effeet? 

A likely important elue is eonservation of magnetie helieity, 

^{A-B) + V-{c(pB + cExA) = -2cE-B. (11) 

at 

This equation, derived direetly from Maxwell’s equations (where cp and A are the 
standard eleetromagnetie sealar and veetor potentials), shows that the volume integral 
of the “magnetie helieity density” A - B is eonserved if E ■ B = 0, whieh is the ease 
whenever the resistivity is zero, so the generalized Ohm’s law reduces toE-|-vx5/c = 0, 
regardless of the specific form of the velocity field v. In partieular, the stable MHD 
equilibria found by Braithwaite & Spruit |0|] most likely minimize energy at a given 
magnetie helieity, and the Hall effeet also eonserves this quantity. 

Sinee, dimensionally, the magnetie helieity density is ~ B^L (where B and L are a 
eharaeteristie magnetie field strength and length seale), whereas the magnetie energy 
density is ~ the helieity tends to reside on large seales, making it mueh more 
diffieult to dissipate than energy. This eould make strongly helieal eonfigurations stable 
and prevent them from deeaying under any of the proeesses eonsidered (aside from 
the extremely slow resistive diffusion). We are currently engaged in identifying such 
configurations. 


CONCLUSIONS 

The evolution of the magnetie field in neutron stars is a very ehallenging subjeet, with 
strong observational elues and involving eomplex physios. Perhaps surprisingly (given 
the extreme oonditions of density, gravity, and field strength), it is likely to involve many 
of the same proeesses also emerging in the evolution of magnetie fields in other oontexts, 
from plasma physios laboratories to galaxies, so oontinued oross-feeding of insights 
should be benefioial both to the study of neutron stars and of other systems. 
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